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We search for the decay B° — + p'^p" in a data sample of about 227 million T{AS) BB decays 
collected with the BABAR detector at the PEP-II asymmetric-energy e^e~ collider at SLAC. We 
find no significant signal and set an upper limit of 1.1 x 10~^ at 90% CL on the branching fraction. 
As a result, the uncertainty due to penguin contributions on the CKM unitarity angle a measured 
inB ^ pp decays is 11° at 68% CL. 

PACS numbers; 13.25.Hw, 11.30.Er, 12.15.Hh 



Measurements of CP-violating asymmetries in the 
B^B^ system provide tests of the standard model by over- 
constraining the Cabibbo-Kobayashi-Maskawa (CKM) 
quark-mixing matrix through the measurement of 
the unitarity angles. Measuring the time-dependent CP 
asymmetry in a neutral- i?-meson decay to a CP eigen- 
state dominated by the tree-level amplitude h — > uud 
gives an approximation a^s to the CKM unitarity angle 
a = arg [-V^jV^jV^^V*^]. The correction Aa = a-acff, 
which accounts for the effects of penguin-amplitude con- 
tributions as an additional decay mechanism, can be ex- 
tracted from an isospin analysis of the branching frac- 
tions of the B decays into the full set of isospin-related 
channels j^j. 

Measurements of branching fractions and time- 
dependent CP asymmetries in B ^ tttt, p7r, and pp have 
already provided information on a. Because the branch- 
ing fraction for B^ tt^tt^ is comparable to that for 
_B+ tt+tt" and tt+tt^, the limit on the cor- 

rection is weak: |Aa^7r| < 35° at 90% confidence level 
(CL) 0. (Charge conjugate B decay modes are im- 
plied in this paper.) In contrast, the channel has a 
much smaller branching fraction than the channels with 
charged p's 0, IE Hi U As a consequence, it is possi- 
ble to set a tighter limit on Aopp. This makes the pp 
system particularly effective for measuring a in a model- 
independent way. 

In B ^ pp decays the final state is a superposition 
of CP-odd and CP-even states, and an isospin-triangle 
relation j3| holds for each of the three helicity amplitudes, 
which can be separated through an angular analysis. The 
measuredpolarizations in B^ P^p'^ ^ and B^ 
P^P~ |y>L3| modes indicate that the p's are nearly entirely 
longitudinally polarized. The current best limit on the 
^0 pQpQ branching fraction was obtained by BABAR 
with a sample of 89 million T{AS) BB decays |^] . 

In this Letter we present improved constraints on the 
^0 pOpO bj-ancJiing fraction and the penguin contribu- 
tion to the measurement of the unitarity angle a. These 
results are based on data collected with the BABAR de- 
tector Q at the PEP-II asymmetric-energy e+e^ col- 
lider located at the Stanford Linear Accelerator Cen- 
ter. A sample of 226.6 ± 2.5 million BB pairs, corre- 
sponding to an integrated luminosity of approximately 
205 fb^^, was recorded at the T(4S') resonance with 




FIG. 1: Definition of helicity angles 6i, 82, and (j) for the 
decay — > p^p". The final states are shown in the p" 
rest frames. 

the center-of-mass (cm.) energy y/s = 10.58 GeV. 
We use a sample of 16 fb~^ taken 40 MeV below the 
T{4:S) resonance to study background contributions from 
e+e^ qq {q ~ u, d, s, or c) continuum events. 

To reconstruct B^ p^p^ {tt'^tt~){tt~^tt~) candi- 
dates, we select four charged tracks that are consistent 
with originating from a single vertex near the e+e^ inter- 
action point. Particle identification is provided by mea- 
surements of the energy loss in the silicon vertex tracker 
and the drift chamber and by the Cherenkov angle in an 
internally refiecting ring-imaging Cherenkov detector. 

The angular distribution of the B'^ p^ p^ decay prod- 
ucts can be expressed as a function of the helicity angles 
(01, 6*2, (f), which are defined by the directions of the 
two-body p^ decay axes and the direction opposite the 
B in the p° rest systems, as shown in Fig. ^ Since the 
detector acceptance does not depend on 4>, the resulting 
angular distribution d^V /{T dcosOi dcos02) is 

^ I J (1 - /l) sin^ 01 sin2 02 + /l cos^ 6^ cos" 02 1 , (1) 

where Jl = |^oP/(2^|^Ap) is the longitudinal polariza- 
tion fraction and Aa=„i.o,+i are the helicity amplitudes. 

The identification of signal B candidates is based on 
two kinematic variables: the beam-energy-substituted 
mass, mEs = [(s/2 + Pi • PB)'^/Ef - p|]^/^, where the 
initial total e'^e~ four-momentum (Ei, pi) and the B mo- 
mentum pb are defined in the laboratory frame; and the 
difference between the reconstructed B energy in the cm. 
frame and its known value AE = i?™ — ^/s/2. The signal 
mES and AE resolutions are 2.6 MeV/c^ and 20 MeV, 
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respectively. The selection requirements for niES, ^E, 
the two 7r~'"7r~ invariant masses to 1,2, and the helicity 
angles are the following: 5.24 < toes < 5.29 GeV/c^, 
|A£;| < 85 MeV, 0.55 < mi,2 < 1-00 GeV/c^, and 
I cos 6*1^2 1 < 0.99. The latter requirement removes a re- 
gion with low reconstruction efficiency. 

To reject the dominant continuum background we re- 
quire I cos^tI < 0.8, where 6t is the angle between the 
B-candidate thrust axis and that of the remaining tracks 
and neutral clusters in the event, calculated in the cm. 
frame. Other discriminating variables include the polar 
angles of the B momentum vector and the B-candidate 
thrust axis with respect to the beam axis in the cm. 
frame, and the two Legendre moments Lq and Lp of 
the energy flow around the i?-candidate thrust axis . 
These variables are combined in a neural network, the 
output of which is transformed into a variable £ for 
which the signal and background distributions are ap- 
proximately Gaussian. 

We veto the background mode B'-* D^tt^ 
/i^TT^TT^TT^, where ft,^ refers to a pion or kaon. We re- 
quire the invariant mass of the three-particle combina- 
tion that excludes the highest-momentum track in the 
candidate B rest frame to be inconsistent with being the 
Z?-meson mass {\mhiT-R — tod| > 13 MeV/c^). After ap- 
plication of all selection criteria, iVcand = 35740 events 
are retained, most of which are background events with 
candidates in the sidebands of the variables. On average 
each selected event has 1.05 candidates. When more than 
one candidate is present in the same event, one candidate 
is selected randomly. 

The signal selection efficiency determined from Monte 
Carlo (MC) [ill simulation is 27% or 32% for longi- 
tudinally or transversely polarized events, respectively. 
MC simulation shows that 22% of longitudinally and 
8% of transversely polarized signal events are misrecon- 
structed with one or more tracks not originating from the 
B^ decay. These are mostly due to combinatorial 

background from low-momentum tracks from the other 
B. We treat these as part of the signal. 

Further background separation is achieved by the use 
of multivariate i3-flavor-tagging algorithms trained to 
identify primary leptons, kaons, soft pions and high- 
momentum charged particles from the other B in the 
event [T3 |. The discrimination power arises from the 
difference between the tagging efficiencies for signal and 
background in the five tagging categories ctag. 

We use an unbinned extended maximum likelihood fit 
to extract the i?" p^ p^ event yield. The likelihood 
function is 

where Uj is the number of events for each hypothesis 
j (signal, continuum, and six _B-background classes). 



and Vj{xi) is the corresponding probability density 
function (PDF), evaluated with the variables Xi = 
{ttiesj Ai?, £, mi, TO2, cos 6*1, cos 6*2, Ctag} of the «th event. 

We use MC-simulated events to study the background 
from other B decays. The charmless modes are grouped 
into five classes with similar kinematic and topological 
properties: o^^tt^F; B° p°K*°\ B+ p+p°; a 

combination of B ^ pir and 5° — > p~^ p~] and B decays 
to other charmless modes not included explicitly. One 
additional class accounts for the remaining neutral and 
charged B decays to charm modes. The number of events 
in each class rij is left free in the fit with the exception of 
three classes where rij is fixed either to the expectations 
from independent measurements (78 ±20 events of 5+ 
p^ p^ , and 48±8 events of B ^ pn and B'^ P^p ) or to 
the extrapolation from the flavor-SU(3)-related B-decay 
modes (25 ± 18 events of B° p°K*°). 

Since the correlations among the variables are found to 
be small, we take each Vj as the product of the PDFs for 
the separate variables. Exceptions are the correlation be- 
tween the two helicity angles in signal, and mass-helicity 
correlations in backgrounds and misrecontructed signal, 
taken into account as discussed below. 

We use double-Gaussian functions to parameterize the 
Toes and AE PDFs for signal, and a relativistic P-wave 
Breit-Wigner (BW) convoluted with a Gaussian resolu- 
tion function for the resonance masses. The angular dis- 
tribution for signal, expressed as a function of the longi- 
tudinal polarization in Eq. is multiplied by a detector 
acceptance function ^(cos 0i, cos ^2), obtained with MC 
simulation. The distributions of misreconstructed signal 
events are parameterized with empirical shapes in a fash- 
ion similar to that used for B background, as described 
below. The £ variable is described by two asymmetric 
Gaussian functions with different parameters for signal 
and background distributions. 

The Toes distribution of the continuum bac kgro und is 
described with the ARGUS parameterization (lj|. The 
Ai? and resonance mass toi.2 PDFs are parameterized 
with low-degree polynomials. The parameterization of 
the TOi and TO2 distributions includes a BW resonant 
component to account for the real p° resonances in the 
continuum background, which are assumed to be unpo- 
larized and thus to have a flat distribution in cos ,2 • The 
cos Oi^2 distribution of the continuum background exclud- 
ing the real resonances is parameterized with a second- 
degree polynomial and an exponential function to allow 
for the increased fraction of combinatorial tt+tt^ candi- 
dates with low momentum pions near | cos ^1^2! = 1- This 
parameterization depends on the p candidate's mass. 

The PDFs for exclusive non-signal B decay modes are 
generally modeled with empirical non-parametric distri- 
butions 0]. However, analytical distributions are used 
for the variables that have distributions identical to those 
for signal, such as toes when all four tracks come from 
the same B, or ■n^'K~ invariant mass TOi^2 when both 
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tracks come from a meson. The two candidates of 
some exclusive non-signal modes can have very different 
mass and helicity distributions. This occurs when one of 
the two p^ candidates is real (e.g., p'^ p^ , p'^K*^) or when 
one of the two p° candidates contains a high-momentum 
pion (ai7r). In such cases, we use a four- variable corre- 
lated mass-helicity PDF. 

The signal and ^-background PDF parameters are ex- 
tracted from MC simulation while the continuum back- 
ground PDF parameters are obtained from data in tties 
and Ai? sidebands. The MC parameters of mEs, 
and £ are adjusted by comparing data and MC in cal- 
ibration channels with similar kinematics and topology, 
such as D~7r+ with D~ K^it~it~ . Finally, the 

B-flavor tagging PDFs for all categories are the discrete 
Ctag distributions of tagging efficiencies. Large samples 
of fully reconstructed i?-meson decays are used to obtain 
the i3-tagging efficiencies for signal B decays and to study 
systematic uncertainties in the MC values of B-tagging 
efficiencies for the B backgrounds. 

Table m shows the results of the fit. No significant sig- 
nal yield is observed. We obtain an upper limit by in- 
tegrating the normalized likelihood distribution over the 
positive values of the branching fraction. The value of 
Jl is fixed to 1 in the fit, as this assumption has been 
shown to give the most conservative upper limit and it 
approximates the values obtained in the B pp decays 
dominated by the tree-level amplitude. The statistical 
significance is taken as the square root of the change in 
— 21n£ when the number of signal events is constrained 
to zero in the likelihood fit. In Fig. [21 we show the pro- 
jections of the fit results onto toes and AE. 

Systematic errors in the fit originate from uncertain- 
ties in the PDF parameterizations, which arise from the 
limited number of events in the sideband data and sig- 
nal control samples. The PDF parameters are varied by 
their respective uncertainties to derive the corresponding 
systematic errors (6.0 events). The event yields from the 
_B-background modes fixed in the fit are varied according 
to the uncertainties in the measured or estimated branch- 
ing fractions. This results in a systematic error on the 



TABLE I: Summary of results: signal yield (risig), selection 
efficiency (Eff), branching fraction (B), branching fraction up- 
per limit (UL) at 90% CL, and significance (including sys- 
tematic uncertainties). The assumption /l = 1 is used. The 
systematic errors are quoted fast. 



Quantity 



Value 



Wsig (events) 
Eff (%) 
B (xlO-^) 
UL (xlO-«) 
Significance (a) 



33I20 ± 12 
27.1 ± 1.3 



0.54+°o ll ± 0.19 
1.1 
1.6 




5.25 5.26 5.27 5.28 



■0.05 0.05 

AE (GeV) 



FIG. 2; Projections of the muitidimensional fit onto mES and 
AE with a requirement on the signal-to-background proba- 
bility ratio 'Psig/'Pbkg with the plotted variable excluded. The 
histogram shows the data and the solid (dashed) line shows 
the full (background only) PDF projection. The projections 
contain 22.5% and 23.9% of signal, and less than 0.5% and 
0.2% of continuum background, respectively. 



signal yield of 5.8 events. We also assign a systematic 
error of 3.0 events to cover a possible fit bias, evaluated 
with MC experiments. 

We estimate the systematic uncertainty due to signal- 
afn^ interference using a simulation study in which the 
decay amplitudes for 5° 
ing to this measurement and those for B'^ 



p^ p° are generated accord- 



cor- 



respond to a branching fraction of 4 x 10^^ (15|]. The 
relative phases between these are modeled with BW am- 
plitudes for all p — !■ TTTT and ai — )■ pi: combinations, with 
additional constants. The values of the constants and 
the a^TT^ CP asymmetries are varied over the allowed 
ranges. We take the rms variation of the average signal 
yield (7.5 events) as a systematic uncertainty. 

Uncertainties in the reconstruction efficiency arise from 
track finding (3%), particle identification (2%), and other 
selection requirements, such as on vertex probability 
(2%), track multiplicity (1%), and thrust angle (1%). 

Our measurement confirms the small value of the 
B^ p^p^ branching fraction with the statistical uncer- 
tainty improved by approximately a factor of two over 
our previous result . Since the tree contribution to the 
— > decay is color-suppressed, the decay rate is 
sensitive to the penguin amplitude. Thus, this mode has 
important implications for constraining the uncertainty 
due to penguin contributions in the measurement of the 
CKM unitarity angle a with B ^ pp decays. 

In the isospin analysis Q, we minimize a that in- 
cludes the measured quantities expressed as the lengths 
of the sides of the isospin triangles. We use the measured 
branching fractions and fractions of longitudinal polar- 
ization of the B+ p+p^ 110] and B° p+p- HQ 



and C+- 



ob- 



decays, the CP-violating parameters 
tained from the time evolution of the longitudinally po- 
larized P^ P^ decay and the branching fraction 
of 5° 



pOpO fj.Qj-Q ^j^^jg analysis. 



We neglect isospin- 
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-10 , 0, . , 10 
Aapp (deg) 



FIG. 3: Ax^ on Aopp obtained from the isospin analysis 
discussed in the text. The dashed lines at Ax^ ~ 1 and 
Ax^ = 2.7 are taken for the 1 a (68%) and 1.64 a (90%) 
interval estimates. 
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breakin g eff ects, non-resonant, and 1=1 isospin contri- 
butions^. 

With the p'^p'^ measurement we improve the con- 

straint on a due to the penguin contribution and obtain a 
68% (90%) CL limit on Aapp = a - aoff of ±11° (±14°). 
Fig. 01 shows the Ax^ on Aapp. Since the central value 
from Fig. 121 is Aapp = 0, the central value of a obtained 
from the isospin analysis is the same as a^B, which is con- 
strained by the relation sin(2acff) = S^~/{1 - Ct"^)^/^ 
and is measured with the B^ P^P^ decay |7[ to be 
ckeff = (I02ti^(stat)t^(syst))° at 68% CL, where the so- 
lution closest to the CKM best fit central value |0| is 
chosen. 

The error due to the penguin contribution may become 
the dominant uncertainty in the measurement of a using 
B => pp decays. However, if B^ p^p^ decays are ob- 
served, time-dependent and angular analyses will allow 
us to measure the CP parameters S^^ and C^", analo- 
gous to and , resolving ambiguities inherent to 
isospin triangle orientations. 

In summary, we have improved the precision on the 
measurement of the —> p^p^ branching fraction by ap- 
proximately a factor of two. The limit on this branching 
fraction relative to those for B+ P^ and B^ P^P~ 
provides a tight constraint on the penguin uncertainty in 
the determination of the CKM unitarity angle a. The re- 
sults summarized in Table |3 supersede our previous mea- 
surement |4. 
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